Spin- wave dispersion in half-doped La3/2Sri/2 Ni04 
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Recent neutron scattering measurements reveal spin and charge ordering in the half-doped nick- 
elate, La3/2 Sri/2 Ni04. Many of the features of the magnetic excitations have been explained in 
terms of the spin waves of diagonal stripes with weak single-ion anisotropy. However, an optical 
mode dispersing away from the (7r,7r) point was not captured by this theory. We show here that 
this apparent optical mode is a natural consequence of stripe twinning in a diagonal stripe pattern 
with a magnetic coupling structure which is two-fold symmetric, i. e. one possessing the same spatial 
rotational symmetry as the ground state. 
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Strongly correlated electronic systems often exhibit 
some evidence of charge, spin, or orbital order, or some 
combination thereof. The interplay between these de- 
grees of freedom can lead to a wide variety of novel 
phases. The nickelates in particular show a wide range 
of doping in which both charge and spin order coexist 
in stripe patterns. The less-ordered stripe structures in 
some families of cuprates have been widely studied for 
their possible connection to high temperature supercon- 
ductivity. 

Recent experiments on doped Ija,2-x Sr^ Ni04 have 
shown clear evidence of static diagonal charge and spin 
stripe order, [ij, ,2i Spin wave theory has been successful 
at describing much of the behavior in these spin-ordered 
systems. 0, ^ IE IE IE IS We consider here the recent 
experiments by Freeman et al. on the spin dynamics 
of half-doped La3/2 Sri/2Ni04 using inelastic neutron 
scattering. In this material, the spins are in a diag- 
onal stripe phase, where stripes run 45° from the Ni-0 
bond direction, and the charged domain walls are only 2 
lattice constants apart. The charge density modulation 
can either be considered as densely packed stripes, or as 
a checkerboard, for domain walls centered on the Ni sites, 
since in that case the two ways of describing the charge 
pattern are indistinguishable at this filling. However, if 
the domain walls are centered on oxygen sites {i.e. for 
bond-centered stripes), the charge checkerboard descrip- 
tion is not applicable. 

In this paper, we are interested in the extra mag- 
netic mode dispersing away from the antiferromagnetic 
wavevector Qaf = (0.5, 0.5) above 50meV in Fig. 3 of 
Ref. IE One explanation put forth by the authors is that 
diagonal discommensurations in the spin order may be 
able to account for this extra scattering mode. We show 
here that the mode could also be due to asymmetry in the 
spin coupling constants, in which case the "extra mode" 
is really an extension of the acoustic band, made visible 
due to stripe twinning. 

The observed ordering vector is Q — (0.275,0.275), 
which is close to a commensurate stripe value of Q = 
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FIG. 1: (Color online) Spin-charge ordering in a Ni02 square 
lattice. Arrows represent spins on Ni^"'' sites and circles rep- 
resent Ni^^ holes, (a) DS2: Diagonal, site-centered stripes of 
spacing p = 2. Note that in this configuration the charges dis- 
tribution is symetric under 90° rotations, and is equivalent to 
a checkerboard pattern. The spin configuration is symmetric 
under 180° rotations, and breaks the symmetry of the charge 
checkerboard. The straight-line exchange coupling across the 
charge domain wall is Jb, and Jc and are diagonal couplings 
perpendicular and parallel to the charge domain wall, respec- 
tively, (b) DB2: Diagonal, bond-centered stripes of spacing 
p = 2. Here, both the charge and spin patterns are 2-fold 
rotationally symmetric. 



(0.25,0.25). The slight deviation from commensurabil- 
ity is believed to be due to the discommensurations de- 
scribed above. Since we are interested in describing rel- 
atively high energy effects, in what follows we neglect 
the small incommensurability, and consider commensu- 
rate diagonal stripe structures of spacing p = 2. We con- 
sider two patterns in this paper: site-centered stripes as 
shown in Fig.^a), and the corresponding bond-centered 
stripes shown in Fig.^b). 

To model these two systems, we use a suitably 
parametrized Heisenberg model on a square lattice. 
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where the indices i and j run over all sites, and the cou- 
plings Jij are illustrated in Fig. ^ For diagonal, site- 
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FIG. 2: Band structures of diagonal site-centered stripes 
(DS2). (a) Jc = Jd = 0; (b) Jc = Ji, Ja = -0.5Ji,. 

centered stripes of spacing p — 2 (DS2), there is no need 
for nearest-neighbor couphng, and so we set Jq = 0. 
The straight-Hne couphng Jf, across the domain wall is 
antiferromagnetic. The diagonal coupling Jc across the 
domain wall is also antiferromagnetic, but the diagonal 
coupling Jd parallel to the stripes we take to be ferro- 
magnetic, Jd < 0, as explained below. In the diagonal 
bond-centered case (DB2), the nearest neighbor coupling 
Ja > is finite and antiferromagnetic. We also include 
the ferromagnetic coupling Jb < across the domain 
wall. Since we are interested in describing high energy 
effects, we neglect the very weak single-ion anisotropy 
term, which splits the 2-fold degenerate acoustic band at 
low energy, with one mode remaining gapless at the IC 
point Qic — (0.25,0.25), and the other mode developing 
a small gap. 8] 

For the case of site-centered stripes in the absence of 
the diagonal couplings, i.e. for Jc ^ Jd = 0, the spin 
system reduces to two interpenetrating antiferromagnets, 
with two separate Neel vectors but identical Neel order- 
ing temperatures. Any weak diagonal coupling is suffi- 
cient to establish a unique relative direction between the 
two Neel vectors, and the ground state becomes the stripe 
structure shown in Fig. ^a). The number of reciprocal 
lattice vectors is also decreased by a factor of 2 in the 
presence of the diagonal couplings Jc 7^ or 7^ 0, as 
can be seen in the bandstructure of Fig. [3 In cither case, 
independent of the value of Jc and J^, although the anti- 
ferromagnetic point QxF is a magnetic reciprocal lattice 
vector and therefore must have a spin wave cone dispers- 
ing out of it, there is no net antiferromagnetism in the 
system, so that weight is forbidden at zero frequency at 
Qaf- The cone that emanates out of Qaf gains finite 
weight as energy is increased, but remains faint at low 
energies. 

Another key feature of nonzero couplings Jc and Jd 
for site-centered stripes is the symmetry of the spin 
wave structure, as shown in Fig. |21 In the limit where 
Jc = Jd = 0, the spin wave dispersion is symmetric un- 
der 90° rotations, as shown in Fig.|2Ia). However, when 



u) 4 




0.25 0.5 0.75 1 0.25 0.5 0.75 1 

K/2Tt fci/27r 



FIG. 3; DS2: Diagonal, site-centered stripes of spacing p = 
2. The plots are for twinned stripes, summing contributions 
along {kx,kx) and {kx,-kx). (a)Jc = Jt and Jd = — 0.5J(,; 
(b)Jc = 2Jj, and Jd = — O.SJ;,. 

either Jc or Jd or both are nonzero, the symmetry is bro- 
ken, and the spin wave structure is now only symmetric 
under 180° rotations, as shown in Fig. [21[b). This means 
that for any nonzero Jc or Jd, the spin wave velocity of 
the acoustic mode dispersing out of Qaf = (0.5,0.5) is 
different parallel and perpendicular to the stripe direc- 
tion. In the presence of stripe twins, the two velocities 
will appear as two branches in plots oi uj vs. k, as shown 
in Fig. 13 

The magnon dispersion from Eq. ^ can be solved an- 
alytically, and for DS2 we find that 

to{k^,ky)^2^/A^-B^, (2) 

where 

A = 2Jb + Jc-Jd + Jd C0s(fc:r - ky) 

B ~ Jb cos(2A;a;) -|- Jb cos{2ky) + Jc cos{kx + ky). (3) 

There are two different spin wave velocities for the cones 
emanating from the IC peak Qic = (0.25, 0.25) and 
symmetry-related points, one corresponding to spin wave 
velocities perpendicular to the direction of the domain 
walls {i.e. perpendicular to the stripes), and the other 
parallel to the direction of the domain walls, 

= A{2Jb + Jc) 
v\\ = Wi2Jb + Jc)i2Jb-Jd) ■ (4) 

In Fig. 13 we show the expected dispersions and scat- 
tering intensities for DS2. Plots are shown for twinned 
stripes, summing the contributions parallel and perpen- 
dicular to the stripe direction, i.e. along the {kx,—kx) 
and {kx,kx) directions, respectively. Because Jc and Jd 
are nonzero, there is an apparent optical mode. Note 
that it is not a true optical mode, since in this config- 
uration there are only two spins per unit cell, leading 
to only one (twofold degenerate) band. As in the bond- 
centered case, weight is forbidden at low energy at the 
antiferromagnetic point Qaf, so that the spin- wave cone 
emanating from this magnetic reciprocal lattice vector 
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FIG. 4: (Color online) Constant energy cuts with windows of 
O.lJbS for twinned diagonal, site-centered stripes of spacing 
p = 2 a.t Jc — 2Ji, and Jd = — 0.5Jft. The energy E is in units 
of JbS. 



has vanishing weight as w ^ 0. In Fig.l^Ia), we have set 
Jc = Jb and Jd — —0.5 Jb- In Fig. Ol^b), we use Jc — 2Jb 
with Jd = — 0.5 Jfc. Notice that in panel (a) of the figure, 
the apparent optical mode is flat. 

We have chosen the coupling constants with the 
following in mind: The "acoustic" branch peaks at 
w(37r/4,37r/4) = 4Jf, + 2Jc. The apparent "optical" 
mode peaks at a;(7r/2, 37r/2) = 4-^/ (2 Jf, — Jd){Jc — Jd)- 
The data indicate that the apparent optical mode is 
higher in energy than the top of the "acoustic" part: 
uj{tt/2,3tt/2) > w(37r/4, 37r/4), which implies that 



Jd < (l/2)(2 + Jc- V2v/4 + J2) 



(5) 



when Jfc = 1. However, the extra mode above 50meV is 
not too high in energy, so parameters need to be chosen 
so as to satisfy this constraint, but remain close to the 
equality. 

We also require the apparent "optical" branch to be 
concave, since there is no evidence of a dip in the extra 
mode. This requirement gives 



.._4(J,-2)(Jc-2-2J,)_ 
dkl ^ V-(-^c- Jd)(Jrf-2) - 
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FIG. 5; DB2: Diagonal, bond-centered stripes of spacing p = 
2. (See Fig.^b).) The plots are for twinned stripes, summing 
contributions along {kx,kx) and {kx,—kx)- (a) Very weak 
coupling across the domain wall, with Jb — — O.lJa. (b) Weak 
coupling across the domain wall, with Jt — — 0.5 Ja- 



resulting in the constraint that 
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As long as the second constraint, Eqn.[3 is satisfied in 
the range < Jc < 2 J;,, then the first constraint, Eqn.|3 
is also satisfied. To describe the data, then, we find that 
we need a sizeable Jc, on the order of Jb. This makes Jc 
significantly larger than that reported at the lower doping 
X = 1/3, where diagonal site-centered stripes of spacing 
p — 3 (DS3) were used to explain the data successfully. (3| 
Although the fits in Ref. l2 were good for Jc = and 
were not significantly improved by letting Jc increase to 
Jc ~ 0.5Jb, the data were still well described using a 
nonzero Jc- We find that the data at doping a; = 1/3 are 
also well described by taking Jc to be as large as Jb or 
even 2 Jf, as in our Fig. |31 

We also find that we need Jd < 0, i.e. the diagonal 
coupling parallel to the stripes needs to be ferromagnetic, 
in order to describe the data. If one considers the diag- 
onal spin couplings Jc and Jd to be derived from, e.g., 
a perturbative treatment of a single-band Hubbard or 
three-band Emery model on a square Ni-0 lattice, we 
expect Jc = Jd. Given that the spin ground state breaks 
the 4-fold rotational symmetry of the square lattice to 
only 2-fold rotational symmetry, any finite spin-lattice 
coupling results in the two diagonal directions being in- 
equivalent, and leads to Jc ^ Jd. For weak rotational 
symmetry breaking of the square lattice in the perturba- 
tive regime, one expects Jc = Jo + with Jd = Jo — ^ 
where e is small compared to Jo, so that the anisotropy 
between the two diagonal coupling directions is weak as 
well. (This preserves the spin ground state of Fig.^a).) 
We find, however, that this regime of the coupling con- 
stants leads to the apparent optical mode being too low 
in energy to capture the data. This may indicate that 
the materials are far from the perturbative limit of the 
single-band Hubbard or three-band Emery model. 

Fig.^shows constant energy cuts for DS2 correspond- 
ing to the parameters in Fig. |3Ib). An important fea- 
ture of this configuration is that although the antiferro- 
magnetic point Qaf = (0.5, 0.5) is a magnetic reciprocal 
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lattice vector, zero-frequency weight is forbidden there 
by symmetry, since the stripes have no net Neel vector 
at Qaf- Combined with the fact that there is no opti- 
cal band, the DS2 configm'ation can never have spectral 
weight at the antiferromagnetic point Qaf, even at finite 
frequency. Notice that as energy is increased in Fig. ^ 
a faint spin wave cone emerges from Qaf in a light ring 
of scattering, but none of the constant energy plots show 
any weight at Qaf- This is consistent with the constant 
energy plots for La3/2Sri/2 Ni04 shown in Ref. By 
contrast, the corresponding bond-centered configuration 
(DB2) shown in Fig.^b) has an optical band which dis- 
plays a saddlepoint at Qaf, and rather large scattering 
intensity at finite frequency at Qaf as a result. (See 
Fig. El of this paper, as well as Fig. 8 of our previous 
paper, Ref.|J.) 

While diagonal, site-centered stripes of spacing p — 2 
(DS2) are able to account for the behavior of the appar- 
ent optical mode observed to disperse away from Qaf in 
Fig. 3 of Freeman et al. |S| , there are two other high en- 
ergy features which this model has not captured. One 
is the asymmetry in intensity observed above 30meV 
in the spin wave cones emanating from the main IC 
peaks. The other is a mode in the 31 — 39meV range 
propagating away from (h, k) structural reciprocal lat- 
tice points. These (as well as the apparent optical mode) 
have been attributed to discommensurations in the mag- 
netic order. 

In Fig. El we show the expected dispersions and in- 
tensities for the diagonal, bond-centered stripes (DB2) 
shown in Fig. Qlb). This configuration has a true op- 
tical mode. Fig. Ela) shows weak coupling across the 
charged domain walls, with Jb — — O.lJa, and Fig. |3{b) 
has somewhat stronger coupling across the domain walls, 
with Jf, = — O.SJq. Note that in the bond-centered case, 
couplings across the domain walls are ferromagnetic. Re- 
sults are shown for twinned stripe patterns, summing the 
contribution parallel and perpendicular to the stripe di- 
rection, i.e. along the {kx,^kx) and {kx,kx) directions, 
respectively. Although there is a reciprocal lattice vec- 
tor at Qaf, weight is forbidden there at zero frequency, 
since the Neel vector switches sign across the domain 
walls. We have reported the analytic form of the spin 
wave dispersion in this case in a previous publication. ,2,] 

Because this spin configuration is only 180° symmetric, 
the spin-wave velocity emanating from Qaf is different 
parallel and perpendicular to the stripe direction. How- 
ever, the branch emanating from Qaf in the direction 
parallel to stripes has so little weight as to be effectively 
invisible in the plots. This configuration displays a true 
optical mode because there are four spins in the unit cell. 
The optical mode has a saddlepoint at Qaf and finite en- 
ergy, with increased weight at the saddlepoint. For weak 
coupling across the domain walls (| Jb| < \Ja\), the op- 
tical mode always displays significant weight at Qaf at 
finite frequency. This is not supported by the dataQ, 



which at the energies measured display no scattering at 
Qaf and finite frequency. This likely indicates that the 
domain walls are not bond-centered in this material, but 
are probably site-centered. 

In conclusion, we have used linear spin-wave theory 
to describe the magnetic excitations recently observed 
in neutron scattering^ on La3/2Sri/2 Ni04. Many fea- 
tures of the data were captured in the spin wave analysis 
of Ref. 0. Other features, including an apparent optical 
mode dispersing away from Qaf above 50meV, were at- 
tributed to discommensurations in the spin order. We 
show here that the apparent optical mode may also be 
captured in linear spin wave theory by using a spin cou- 
pling configuration that preserves the symmetry of the 
spin ground state. Namely, we have shown that diagonal, 
site-centered stripes of spacing p — 2 capture this mode 
when the pattern of couplings is 2-fold symmetric. This is 
because the 2-fold symmetric coupling pattern gives rise 
to two different spin-wave velocities {i.e. W|| ^ v±) ema- 
nating from the antiferromagnetic point Qaf = (0.5, 0.5). 
For twinned samples, the two velocities are simultane- 
ously visible, and the higher velocity mode W|| parallel to 
the stripes is responsible for the "extra" scattering above 
50meV. Furthermore, this configuration is forbidden to 
display scattering at the antiferromagnetic point Qaf, 
whereas bond-centered stripes have a true optical mode 
with significant weight at Qaf, which is not supported 
by the data. We therefore conclude that the magnetic 
excitations observed in Ref. |^ are consistent with site- 
centered stripes, but not with bond-centered stripes.. 

It is a pleasure to thank D. K. Campbell and 
A. T. Boothroyd for helpful discussions. This work was 
supported by Boston University (DXY), and by the Pur- 
due Research Foundation (EWC). EWC is a Cottrell 
Scholar of Research Corporation. 



[1] A. T. Boothroyd, D. Prabhakaran, P. G. Freeman, S. J. S. 

Lister, M. Enderle, A. Hiess, and J. Kulda, Phys. Rev. B 

67, 100407(R) (2003). 
[2] H. Woo, A. T. Boothroyd, K. Nakajima, T. G. Perring, C. 

Frost, P. G. Freeman, D. Prabhakaran, K. Yamada, and 

J. M. Tranquada, Phys. Rev. B 72, 064437 (2005). 
[3] E. W. Carlson, D. X. Yao, and D. K. Campbell, Phys. 

Rev. B 70, 064505 (2004). 
[4] D. X. Yao, E. W. Carlson, and D. K. Campbell, Phys. 

Rev. B 73, 224525 (2006). 
[5] J. M. Tranquada, P. M. Gehring, G. Shirane, S. Shamoto, 

and M. Sato, Phys. Rev. B 46, 5561 (1992). 
[6] P. Bourges, Y. Sidis, M. Braden, K. Nakajima, and 

J. M. Tranquada, Phys. Rev. Lett. 90, 147202 (2003). 
[7] F. Kr^iger and S. Scheidl, Phys. Rev. B 67, 134512 (2003). 
[8] P. G. Freeman, A. T. Boothroyd, D. Prabhakaran, C. 

Frost, M. Enderle, and A. Hiess, Phys. Rev. B 71, 174412 

(2005). 



